Cosmic coincidence problem and variable constants of 

physics 

q Mubasher Jamil* Farook Rahaman} and Mehedi K^alam* 

O 

* Center for Advanced Mathematics and Physics, National University of Sciences and Technology, 
D Peshawar Road, Rawalpindi - 46000, Pakistan 

^Department of Mathematics, Jadavpur University, Kolkata - 700032, India 
"'"Department of Physics, Netaji Nagar College for Women, Regent Estate, Kolkata - 700092, India 

February 28, 2009 

O" 1 



Abstract 



CO 
> 

r-H The standard model of cosmology is investigated using time-dependent cosmological con- 

stant A and Newton gravitational constant G. The total energy content is described by the 
modified Chaplygin gas equation of state. It is found that the time-dependent constants cou- 
pled with the modified Chaplygin gas interpolate between the earlier matter to the later dark 

OO energy dominated phase of the universe. We also achieve a convergence of the parameter 

lo — > —1, almost at the present time. Thus our model fairly alleviates the cosmic-coincidence 
problem which demands u = — 1 at the present time. 
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1 Introduction 



The astrophysical observations of supernovae of type la give a convincing evidence of 
a universe undergoing an accelerated phase of its expansion history; rather than going 
to deceleration as was expected theoretically [H El El H]. Similar conclusion has been 
deduced from the observations of anisotropies in cosmic microwave background by the 
WMAP [5j E] which favors a low density, spatially flat (fltot ~ 1) universe filled with 
an exotic vacuum energy containing the maximum portion of the total energy density, 
~ 0.7 [7J. This mysterious dark energy is represented by a barotropic equation of 
state (EoS) p = ujp, where p and p is the pressure and the energy density of dark energy, 
with to < — 1. If uj — — 1, it is called 'cosmological constant', — 1 < c<j < — 1/3 is dubbed 
'quintessence', while it is called 'phantom energy' if uj < —1 [HIE]- All these candidates 
have some fundamental problems: The former one presents discrepancy of 120 orders 
of magnitude between theoretical and empirical results [10] , while quintessence requires 
fine tuning of cosmological parameters for a suitable choice of the potential function |llj . 
Also the phantom energy yields eccentric predictions like 'big rip' and ripping apart of 
gravitationally bound objects [T2l ITS"] . Moreover, the variations in uj suggest that there is 
no consensus on the actual EoS of dark energy and one can only deal with the upper and 
lower bounds on uj [14|. 

Although the dark energy is generically considered to be a perfect fluid yet it need not 
to be perfect if it experiences perturbations [15]. Moreover, the dark energy also may not 
be completely 'dark' especially if it gets coupled with matter and energy exchange takes 
place; consequently the matter evolution becomes modified [T6J |TT] . It is now obvious that 
this sudden transition to acceleration from the earlier deceleration phase is rather recent 
with the corresponding redshift z < 1 [18]. It has been proposed that the standard model 
of cosmology may not be sufficient to explain this exotic phenomenon and hence significant 
modifications are proposed like a modified Friedmann equation H 2 ~ g{p), where g is an 
arbitrary function of p and H is the Hubble parameter [TjJ] and adding a Cardassian 
term in the Friedmann equation [20] • Other possible explanations proposed are dark 
energy arising from tachyonic matter [2T] [22] . van der Waals fluid [23], geometric dark 
energy [21], Dvali-Gabadadze-Porrati (DGP) gravity model [25] and Randall-Sundrum 
brane world model [26J are most prominent. Although the nature of dark energy is not 
clear but its thermodynamical properties suggest a universe filled with it becomes hotter 
with time [27] . 

One of the most interesting problems in the present cosmology is the cosmic coincidence 
(or cosmic conundrum) problem which naively asks: why the energy densities of matter 
and dark energy are of the same order or the corresponding dimensionless ratio is closer to 
1, at current time [2EJ I2HJ EO] • This problem can be posed in terms of the EoS parameter uj: 
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if the parameter u > —1 in the past and to < — 1 in the future then why we are observing 
u — — 1 at present time. In recent years, this problem is addressed using a notion of 
interacting dark energy model in which both the interacting components i.e. dark energy 
and matter exchange energy to keep the density ratio close to 1. This model has some 
intrinsic problems that are still unresolved: the microphysics of energy transfer is not 
exactly understood i.e. the particles that can mediate the interaction are not pointed out. 
Moreover, the coupling function (or the decay rate) for the required interaction is chosen 
quite arbitrarily j3TJ [32} [33j El] and also the coupling constant involved is not yet properly 
constrained theoretically or observationally [35]. In this paper, we address this problem 
using a simplified approach by considering the constants of physics to evolve over cosmic 
time. We here take three ansatz for scale factor and analyze the behavior of parameter 
oj. Curiously, the parameter uj(z) evolves from positive to negative values and finally 
converges to —1. This result turns out to be consistent with the observations. Hence our 
model fairly addresses the cosmic coincidence problem and practically alleviates it. 

Due to multitude of uncertainties in the determination of uo observationally and other 
intrinsic theoretical problems (as discussed above) with it, we here proceed with an EoS 
commonly called the Chaplygin gas (CG) represented by [SU I3T] 



where X is a constant parameter. The CG effectively explains the evolution of the universe 
from the earlier deceleration (matter dominated era) to the later acceleration phase (dark 
energy dominated) as is manifested in the following equation 



matter unification models [38l [39] . Due to its this effectiveness, several generalizations of 
CG are proposed (see e.g. [201 SH SSI SHI SSI SZJ ) - The Chaplygin gas arises from 
the dynamics of a generalized ci-brane in a (d + 1, 1) spacetime and can be described by a 
complex scalar field whose action can be written as a generalized Born-Infeld action |48j . 

The plan of the paper is as follows: In the second section, we shall present the model 
of our system. In third section, we determine the cosmological parameters for different 
choices of the scale factor parameter. The last section is devoted for the conclusion of our 
paper. 




(2) 
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2 The cosmological model with variable constants 



The Friedmann-Robertson- Walker (FRW) metric which satisfies the cosmological princi- 
ple is specified by 



ds 2 = -dt 2 + a 2 {t) 



dr2 — + r 2 {de 2 + ^ 2 

— kr A 



(3) 



Here a(t) is the scale factor that determines the expansion of the universe. Also the 
parameter k is the curvature parameter determining the spatial geometry of the FRW 
spacetime. It can take three possible values k = +1,0, —1 which correspond to spatially 
closed, flat and open universe respectively or geometrically spherical, Minkowskian and 
hyperbolic spacetime respectively. 

The equations of motion corresponding to FRW metric are 

a A-kG , _ . A . 

a 6 6 

Above H is the Hubble parameter. Note that we have assumed k = in the above 
equations which is favored by the observational data. Here A is called the cosmological 
constant with dimensions of (length) -2 . Note that Eq. (5) shows that accelerated expan- 
sion of the universe a > is possible if the strong energy condition p + 3p > is violated 
and also it is independent of the choice of k. The energy conservation equation for the 
above system is 

p + 3H (p + p) = 0. (6) 

The cherished constants of physics that describe the universe need not to be constant 
but can vary with respect to other parameters. For instance, the cosmological constant 
A(t) = A(t ) + (t — t )A(t ) + ... which is constant at zeroth order approximation but is 
really a time dependent function at higher order approximations. Note that accelerated 
expansion of the universe follows from A > 0. The cosmic history of A shows that it 
was large in the past while it is small at present and will continue to decrease, hence it 
gives a parametrization A oc t a , A oc p 1 and A oc H 2 ED]- A variable cosmological 
constant also arises in theories of higher spatial dimensions like string theory and manifests 
itself as the energy density for the vacuum [51 J and it can also addresses the cosmic 
age problem effectively [52]. Similarly, there is some evidence of a varying Newton's 
constant G: Observations of Hulse- Taylor binary pulsar B1913 + 16 gives a following 
estimate < G/G ~ 2 ±4 x 10~ 12 yr" 1 [53J, helioseismological data gives the bound 
< G/G ~ 1.6x 10~ 12 yr -1 [54J (see Ref [55] for various bounds on G/G from observational 
data). The variability in G results in the emission of gravitational waves. Dimensional 
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analysis also shows that the time dependent parameter A to be decreasing with time t 
[56| [57]. In another approach, it is shown that G can be oscillatory with time [58] • It is 
recently proposed that variable cosmic constants are coupled to each other i.e. variation 
in one leads to changes in others (59]. A variable gravitational constant also explains the 
dark matter problem as well [60]. Also discrepancies in the value of Hubble parameter 
can be removed with the consideration of variable G [61J. Due to these reasons, we shall 
take A and G to be time dependent quantities i.e. A = A(t) and G = G(t). Hence Eqs. 
(4) and (5) yield 

Gp + pG + 3(p + p)GH+^ = 0. (7) 
Using Eqs. (6) and (7), we can write 

„G + A = . (8) 
Taking the ansatz for cosmological constant as 



A=^, (9) 

where (3 and 7 are constant parameters. Note that this is a general ansatz and can reduce 
to Chakraborty and Debnath [63] if 7 = — 1. Using the modified Chaplygin gas (MCG) 
EoS given by [M] 

P = A P - j a , (10) 

where A and B are constant parameters and < a < 1. Thermodynamical analysis of 
MCG show that the values a = 1/4 and B = 1/3 are consistent with the phenomenological 
results [65]. It is also shown that the recent supernovae data favors a > 1 values [661 167] . 
The MCG reduces to generalized Chaplygin gas (GCG) if A = while it gives CG if further 
a = 1. While a barotropic EoS is obtained if B — 0. Thus Eq. (10) is a combination of 
a barotropic and GCG EoS. Precisely, the observations of cosmic microwave background 
gives the constraint —0.35 < A < 0.025 at 95% confidence level [68J. Analysis of various 
cosmological models show that models based on Chaplygin gas best fit with supernova 
data |69j. Using Eq. (10) in (6), we get the density evolution of MCG as 

p= (x + Ca^) 1 ^, (11) 

where x = and y = 3(1 + a)(l + A). Making use of Eqs. (9) and (11) in (8), the 
parameter G is determined to be 



l7r(a + 7 + 2) 



[x 



xa y \ s .. I . . . . .r<r< 



+ ca-yy 5 1 + ^- 2^1 [6,6,1 + 5, 



C J ' x v c 



(12) 
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where 

S= a + -> + 2 . (13) 
1 + a v ' 

The parameter A can alternatively be written as 

A = 3(3{x + Ca~ y )^ . (14) 

3 Determination of cosmological parameters 

To analyze the behavior of the above cosmological parameters, we will consider three 
cases: 

1. a ~ a T n , 

2. a ~ (uT-u)^, 



i 



3. a~[e- DPT -l]-i> 

Here Z?, P, q, u, v, n and a Q are constant parameters. Also T — t/t is the dimensionless 
time parameter with t a is the current age of the universe. 



3.1 Power law form of scale factor 

We consider power law form of the scale factor 

a ~ a Q T n , (15) 

where a Q and n are arbitrary constants. For this choice, it is possible to get the accelerated 
expansion of the universe if n > 1. Now, all the physical parameters will take the following 
forms as: 

- 81(0 + 7 + 2)'^ (a T")») \ + C } 

x l F 1 (^,l + 4-^)], (17, 
/ c \ ^ 

A ~ 3/3 x + 7 - r , (18) 
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p ~ A 



C s 
{a T n )y, 



i - 

l+a 



5 



(x + c "l 1+a 



(19) 



The cosmological parameters obtained in this section are plotted in figures 1 to 4 against 
parameter T. 
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Figure 1: The energy density p (see Eq. 16) is plotted against time T. Model parameters are 
fixed at A = B = C = n = 2, a Q = l, 7 = -0.5 and = 1. 
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Figure 2: The parameterized Newton's constant G (see Eq. 17) is plotted against time parameter 
T for different choices of a. Model parameters are fixed as taken in Fig.l 
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Figure 3: The cosmological constant A (see Eq. 18) is plotted against time parameter T. Model 
parameters are fixed as taken in Fig.l 
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a = 1 a = 0.3 a = 0.6 



Figure 4: The EoS parameter uj = p/p (see Eq. 16 and 19) is plotted against the time parameter 
T. Model parameters are fixed as taken in Fig.l 
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3.2 Negative constant deceleration parameter 



In this case, we consider constant deceleration parameter model denned by 

a 

constant = q. 



aH 2 



(20) 



Here the constant is taken to be negative i.e. it is an accelerating model of the universe 
[70] . The solution of equation (20) is 



a ~ (uT - v) 



(21) 



where, u and v are integration constants. This equation implies, the condition of expansion 
is 1 + q > 0. Now, all the physical parameters will take the following forms as: 



C 



x 



G 



(uT - v) P+s 
[ x4 



(22) 



-3/3 7 C 
8n(a + 7 + 2 



x 2 /' i [ (). ). 1 — 5, -x 
C 



c 

(uT - u)Ti+5J 

(uT - v) t 1 ^) 
C 



1 + x 



. y \ 5 

(uT - v) 
C 



A ~ 3/3 | A : 

p ~ A | x + 



(uT - v) 

C 
(uT - 



i 

l+a 



5 



X + 



1 

l+a 



(uT-v) (1+9) 



(23) 
(24) 
(25) 



The cosmological parameters obtained in this section are plotted in figures 5 to 8 against 
parameter T. 
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3.3 Specific form of deceleration parameter 



In this case, we consider specific form deceleration parameter model defined by [7T] 

a Pa p 



1H 2 



l + a 1 



(26) 



Here P is is a constant. This choice of deceleration parameter provides an early deceler- 
ation and late time acceleration of the universe. The solution of equation (26) is 



-DPT 



11"?. 



(27) 



where D is an integration constant. For negative values of P, we always get accelerated 
expansion of the universe. Now, all the physical parameters will take the following forms 
as: 



x + 



C 



i 

l+a 



-j- DPT 



1] 



(28) 



G 



[[x + 



C 



-5 



-DPT 



A 



V 



87r(a + 7 + 2) 
x 2 Fi (s,6,l + 6, -x! 

C 



-DPT 



-DPT 



1] 
1]- 



1 + x- 



p \ 



c 



c 



3/3 [x + 



A\x + 



-DPT 



"I] 



c 



1 

l + a 



B 



-DPT 



I]" 



X 



C 



-1]-p 



l n a ■ 

l+a 



(29) 
(30) 
(31) 



The cosmological parameters obtained in this section are plotted in figures 9 to 12 against 
parameter T. 
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Figure 5: The density parameter p (see Eq. 22) is plotted against time T. Model parameters 
are fixed at A = B = C = 0.08, u = l,v = 2, q = j = -0.2 and (3 = 1. 
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Figure 6: The Newton's constant G (see Eq. 23) is plotted against time parameter T for different 
choices of a. Other model parameters are fixed as in Fig. 5 
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Figure 7: The cosmological constant A (see Eq. 24) is plotted against time parameter T. Other 
model parameters are fixed as in Fig. 5 
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Figure 8: The parameter p/p = u (see Eq. 22 and 25) is plotted against the time parameter T. 
Other model parameters are fixed as in Fig. 5 
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4 Conclusion and discussion 

In this section, we discuss the results of our paper. All the cosmological parameters with 
the exception of u are plotted in figures in 1 to 12 in logarithmic scale against dimen- 
sionless time parameter T. The parameters in section 3.1 are shown in figures 1 to 4. 
The cosmological energy density decreases with time and then remains constant in far 
future. The cosmological constant was large in the past which resulted in inflation while 
now it is small to produce current accelerated expansion. The Newton's gravitational 
constant steadily increased which caused structures to form. Also the dimensionless pa- 
rameter uj varies from the positive to negative values and converging to —1 at current 
time t = t D (~ 1.2H ) or T ~ 1 showing that uj is inherently evolving over cosmic history, 
with T = corresponds to the big bang epoch (units are chosen to be meter, kilogram, 
sec). A similar behavior is obtained for parameters of section 3.2 and 3.3 shown in figures 

5 to 8 and 9 to 12, respectively. In figures (13), (14) and (15), we have plotted the same 
parameters against redshift z. 

The problem attempted in this paper can also be looked in the context of bulk viscous 
cosmology. The anisotropic stresses can be important at large scale and hence they should 
be incorporated in the MCG equation of state (see [72] for the basic formalism). It would 
also be interesting to extend our model using the modified f(R) gravity theory as well 

m. 
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Figure 10: The Newton's constant G (see Eq. 29) is plotted against time parameter T for 
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Figure 13: The energy density p (see Eq. 16) is plotted against redshift parameter z. Model 
parameters are fixed as in Fig. 13 
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Figure 14: The Newton's constant G (see Eq. 17) is plotted against redshift parameter z. The 
model parameters are fixed at A = B = C = 2, 7 = —0.2 and (3=1. 
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Figure 15: The cosmological constant A (see Eq. 18) is plotted against redshift parameter z. 
Model parameters are fixed as in Fig. 13 
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